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ABSTRACT 

Several high-frequency peaked BL Lac objects such as Mrk 501 are strong TeV emitters. However, 
a significant fraction of the TeV gamma rays emitted are hkely to be absorbed in interactions with 
the diffuse IR background, yielding electron-positron pairs. Hence, the observed TeV spectrum must 
be steeper than the intrinsic one. Using the recently derived intrinsic 7-ray spectrum of Mrk 501 
during its 1997 high state, we study the inverse-Conipton scattering of cosmic microwave photons by the 
resulting electron-positron pairs, which implies the existence of a hitherto undiscovered GeV emission. 
The typical duration of the GeV emission is determined by the flaring activity time and the energy- 
dependent magnetic deflection time. We numerically calculate the scattered photon spectrum for different 
intergalactic magnetic field (IGMF) strengths, and find a spectral turnover and flare duration at GeV 
energies which are dependent on the field strength. We also estimate the scattered photon flux in the 
quiescent state of Mrk 501. The GeV flux levels predicted are consistent with existing EGRET upper 
limits, and should be detectable above the synchrotron - self Compton (SSC) component with the 
Gamma-Ray Large Area Space Telescope [GLAST) for IGMFs < 10""'^^ G, as expected in voids. Such 
detections would provide constraints on the strength of weak IGMFs. 

Subject headings: BL Lacertae objects: general — BL Lacertae objects: individual (Markarian 501) — 
diffuse radiation — gamma rays: theory — magnetic fields 



1. INTRODUCTION 

Blazars including high-frequency peaked BL Lac objects 
, (HBLs) are the most extreme and powerful sources among 

■ active galactic nuclei. The standard blazar model con- 
' sists of a supermassive black hole ejecting twin relativis- 
, tic jets, one of which is close to the line of sight. Sev- 

■ eral HBLs such as Mrk 501, Mrk 421, PKS 2155-304, lES 
; 2344-^514, H1426-h428 and lES1959-f 650 are of particu- 
. lar interest because they emit TeV photons (see Catanese 
' & Weekes 1999; Horns et al. 2002). The detection and 
, study of such photons can provide new insights on the en- 
ergetics and physical conditions in the emission regions of 
such blazars (Katarzyhski, Sol & Kus 2001; Kino, Taka- 
hara & Kusunose 2002). Also, constraints on the spec- 
tral energy distribution of the intergalactic infrared back- 
ground can be inferred from the observations on TeV pho- 
tons (for a review see Hauser & Dwek 2001). Stecker, De 
Jager & Salamon (1992) have emphasized that the high- 
energy gamma photon spectra from these blazars will be 
modified by strongly redshift-dependent absorption effects 
due to interactions of such photons with the intergalactic 
infrared- UV background, and indicated that the intrinsic 
spectrum of an observed TeV blazar can be derived by 
evaluating the optical depth to TeV photons. Such cal- 
culations were made for Mrk 501 during the 1997 flaring 
activity, leading to an inferred intrinsic high-energy spec- 
trum with a broad, flat peak that is much higher than the 
observed one in the ~ 5 — 10 TeV range (Konopelko et al. 
1999; De Jager & Stecker 2002, hereafter DS). The physical 
reason for this difference is that a significant fraction of the 
original high-energy gamma rays have been absorbed in 77 
interactions with photons of the intergalactic infrared- UV 



background, leading to electron/positron pairs. 

The purpose of this Letter is to suggest that inverse 
Compton (IC) scattering of the resulting electron/positron 
pairs against cosmic microwave background (CMB) pho- 
tons may produce a new GeV emission component in TeV 
blazars. For gamma-ray bursts, such Compton scatter- 
ing leads to an observable, delayed MeV-GeV emission 
component if the intergalactic magnetic fields (IGMFs) 
are very weak (Plaga 1995; Cheng & Cheng 1996; Dai 
& Lu 2002). A similar phenomenon is also expected from 
gamma-ray burst proton interactions with the CMB (Wax- 
man & Coppi 1996). Here we discuss the well-studied 
blazar Mrk 501, both because the high-energy spectrum 
up to 20 TeV of strong flares of this blazar in 1997 has 
been observed by the HEGRA air Cerenkov telescope sys- 
tem (Aharonian et al. 1999, 2001), and because the in- 
trinsic spectrum of Mrk 501 over two decades of energy 
has been derived based on the consistency between the 
Whipple telescope and HEGRA spectra (DS). 

The strength of IGMFs has not been determined so far. 
Faraday rotation measures imply an upper limit of ^ 10~^ 
G for a field with 1 Mpc correlation length (Kronberg 1994 
for a review). Other methods were proposed to probe 
fields in the range 10"^° G to lO^^o q (Lgg^ Qlinto & Sigl 
1995; Plaga 1995). To interpret the observed /xG mag- 
netic fields in galaxies and X-ray clusters, the seed fields 
required in dynamo theories could be as low as 10~^° G 
(Kulsrud et al. 1997; Kulsrud 1999). Furlanetto & Loeb 
(2001) argued that quasar outflows may pollute the inter- 
galactic medium, but the possible weak IGMFs in voids 
may remain uncontaminated. Theoretical calculations of 
primordial magnetic fields show that these fields could be 



2 



of order 10~^° G or even as low as 10^^® G, generated 
during the cosmological QCD or electroweak phase tran- 
sition respectively (Sigl, Olinto & Jedamzik 1997). In this 
Letter we propose that by observing a hitherto undiscov- 
ered GcV emission component from flares of TcV blazars 
such as Mrk 501, one may be able to obtain important 
information or constraints on the poorly known IGMFs. 

2. PROPERTIES OF EXTERNAL IC EMISSION 

We consider a strong flare, e.g. in Mrk 501, of du- 
ration <var- A fraction of the high-energy 7-rays emit- 
ted during such a flare can be absorbed in the cosmic 
background radiation fields as these 7-rays travel towards 
the observer. The pair production optical depth r®^ de- 
pends strongly on the 7-ray energy (E^) and the redshift 
{z). DS numerically calculated t^^ as a function of the 
photon energy for low redshifts by using the models of 
Malkan & Stecker (2001) of the infrared background radi- 
ation, and extrapolating these models into the optical- UV 
range in terms of recent galaxy count data. DS consid- 
ered both "fast evolution" and "baseline" cases, which 
may be considered to bracket the spectral energy distri- 
bution of the intergalactic infrared-UV background radi- 
ation. If primary photons of energy are absorbed, 
the resulting electron/positron pairs have Lorentz factors 
7e ~ E^/(2meC^) = 10^(£'T,/lTeV), where me is the elec- 
tron mass. The pairs will subsequently Compton scatter 
on the ambient CMB photons. As a result, the initial en- 
ergy of a microwave photon, e, is boosted by IC scattering 
up to an average value ~ 7^6 ~ 0.63(£'-y/l TeV)^ GeV, 
where e = 2.7kT is the mean energy of the CMB photons 
with T ~ 2.73 K and k is the Boltzmann constant. 



to magnetic deflection becomes Ate 
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2.1. The GeV Emission Duration 

Several timescales are involved in the emission pro- 
cess. The first is the well-known angular spreading time. 



i?pair/(27e'c) = 96(7e/106) (nffi/l cm-^) " 



where impair = (0.26c7TniR)-^ ~ 5.8 x 10^4(71^/1 cm'^)-! 
cm is the typical pair-production distance, ctt is the Thom- 
son cross section, and niR is the intergalactic infrared pho- 
ton number density (see Dai & Lu 2002). Therefore, pri- 
mary TeV photons have a typical mean free path of a few 
Mpc, so that although the source may be in a region of 
a high field with ~ 10 ~^ G, these photons may escape to 
much lower field regions. Six TeV blazars detected so far 
have redshifts of ^ 0.03 to ^ 0.1, implying that their lumi- 
nosity distances are much larger than i?pair- Only in such 
cases are our calculations valid. 

The IC cooling timescale (in the local rest frame) of 
relativistic electrons with Lorentz factor of je is tic = 
3mec/(47eO-TMcmb) = 7.3 X 10^3(7e/10'')"^ s, where Ucmb 
is the CMB energy density. Thus the typical flight path, in 
which most of the electron energy is lost, is Aic — ctic = 
2.2 X 102''(7e/10*^)"i cm. This is much less than the dis- 
tance from the source to the observer, implying that en- 
ergy loss of the electron due to IC scattering is local. In the 
absence of any IGMF, the IC cooling timescale in the ob- 
server frame would be Aiic ^ tic/(27e) = 37{je/W^)~^ s. 

In the presence of IGMFs with strength Biq , the elec- 
trons will be deflected. The deflection angle is estimated 
by ^ Aic/iiL = L3 x 10-5(7,/106)-2(Big/10-20G), 
where i?L = lefneC? / {eBio) is the Larmor radius of the 
electrons. The contribution to the emission time due 



6.1 X 10''(7,/10f')--^(BiG/10-20G)2s for 61b < 1. This 
may provide a probe of weak IGMFs and its application 
to the delayed high-energy emission of GRB 940217 yields 
IGMFs which are as weak as Biq ^ 10"^° G. 

Therefore, the duration estimate of the IC emission 
from electron/positron pairs scattering off the CMB is 
Af(7e) = max(Atic, A^A, AfB,ivar)- Taking Ate = ^var, 
which implies a scattered photon energy ~ E'turn = 
0.22(BiG/10~20G)4/5(tvar/lday)-2/5GeV, we find that 
the duration of the expected GeV emission is always given 
by the variability timescale of the TeV gamma-ray flux for 
photon energes larger than i?turn, and it is given by A^b 
for photon energies smaller than £^turn- 

2.2. The Emission Spectrum 

The optical depth (r^^) to high-energy photons was used 
by DS to derive the intrinsic 7-ray spectrum of Mrk 501 
during its 1997 high state, which can be parameterized as 
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where E^ = l(7e/10^) TeV, and the parameters Eb (in 
TeV), K (in 10^^" ergs cm^^ s^^), Fi, r2 and / are shown 
in Table 3 of DS. Please note that equation (1) is a time- 
averaged spectrum of a long-term outburst of about 6 
months in 1997. During this period Mrk 501 has shown 
a number of one-day flares, whose gamma-ray flux values 
appear to be larger than the time-averaged ones by a fac- 
tor of r] with 1 ^ ;^ 3. Thus, for accurate calculations of 
the GeV emission one should use the specific flux values 
for the flare considered. Here we take 7/ = 1 as a conserva- 
tive value. Letting the luminosity distance to the source 
be Dl ~ 5 X 10^^ cm, we have the total electron energy 
spectrum (including the positrons). 
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where C = 57r x 10 °D\K. Therefore, the observed time- 
averaged scattered-photon spectrum is given by (Blumen- 
thal & Gould 1970) 
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where ^ = tvar/At, E~.,j is the externally scattered photon 
energy, t is the time measured in the local rest frame, and 
dN~f^_rJ d,tdE^i^i, expressed by equation (2.48) of Blmnen- 
thaf & Gould' (1970) using the Klein-Nishina (KN) cross- 
section formula, is the spectrum of photons scattered by an 
electron with Lorentz factor of 7e from a segment of the 
CMB photon gas of differential number density n{e)de. 
If the electron energy spectrum is simplified as oc 7^7 '"S 
then a primary analysis similar to Dai & Lu (2002) shows 
that in the Thomson limit the scattered photon spectrum 
oc E-^['+^^/^ for AtB(7e) < ivar and cx E'^^'-^^^^ for 
^tsije) ^ ivar, implying a spectral turnover at scattered 
photon energies of ~ Stum- Figure 1 shows the energy 
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spectra of scattered photons for different assumed fGMFs 
during the f997 flaring activity of Mrk 50f . For an IGMF 
^ 10~^° G, we indeed see a spectral turnover, wliose en- 
ergy is strongly dependent on the field strength. Figure 2 
presents the emission spectra for different flare durations, 
and we find that for a fixed IGMF, the longer the dura- 
tion, the smaller the turnovc;r cncagy becomes and thus 
the emission becomes easier to detect. 

3. COMPARISON WITH THE SSC SPECTRUM 

The observed low-energy hump of Mrk 501 is usu- 
ally thought to be due to synchrotron radiation from 
a relativistic plasma blob. The comoving spectral syn- 
chrotron photon energy density is given by u'{e'^) = 
[2£)|(^(e-y)]/(r^c^^), where (j){€j) is the observed syn- 
chrotron flux at energy of e-y, rb is the comoving radius, 5 
is the Doppler factor, and primes denote comoving quan- 
tities. Here e-y = Se'^, and = (5ctvar- 

In order to calculate the comoving spectral synchrotron 
photon number density n'{e'^), we approximate the flux 



density F{e~^) 



for €^ < Ebr and F(e^ 



for fbr < £7 < Cpk observed during the flares of 1997. 
From the observations by Catancsc ct al. (1997), Plan et 
al. (1998) and Petry et al. (2000), it can be seen that 
ai ~ 0.4, Q!2 — 0.6, ebr — 2 keV and Cpk ~ 100 keV, and 

the synchrotron radiation flux </'(e-y) = ejF{e-y) reaches a 
maximum value (f>p\^ ^ 6 x 10^^° ergs cm^ s^^ at e.y = Cpk- 

3.1. The Intrinsic Optical Depth 

We first derive a generic formula for the optical depth 
to a photon with > ITeV inside the blob. In the 

comoving frame of the blob, this photon can annihilate 
a second softer photon with energy greater than e'^jj = 
d{mec'^)'^ / Ej, yielding an electron/positron pair. The co- 
moving total number density of soft photons that have 
energies greater than ef^„ is calculated by 



2 J, . / ^. \ -"2+1 / ^ 



/— ai 



(4) 

for e!,jj < e\,y/5. Thus, the optical depth is t™^ = 

(11/180)(Tt«'(> ^'anYbi where (11/180)o-t is the averaged 
cross section assuming an F{e~^) oc e'^^ spectrum (Svens- 
son 1987; also see Lithwick & Sari 2001 and Zhang & 
Meszaros 2001). Inserting equation (4) into the above 
equation yields 
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where 



(ll/90)c7T£'i(.^pk/epk)(epk/ebr)"^(ebr/meC^)"^ 
T= . (6) 



For Q!i 



Q2 



1, equation (5) becomes r™ cx i5 ^, which 



is consistent with the previously used optical depth (see 
Aharonian et al. 1999; DS). Inserting the observed val- 
ues of the parameters into equations (5) and (6). we have 
T^^ = 0.015(£;T,/lTeV)0-4(,5/10)-4-8(i^^jlday)'-i. The 
observed highest energy E^ ~ 20 TeV implies that the op- 
tical depth to a photon with such an energy is less than 
unity. This provides a lower limit to the Doppler factor 



5 > 5.4(.B.^./20TcV)l''l2(^var/lday)-l/''■^ i.e., the photon- 
energy cutoff due to pair production in the blob can be 
larger than 20 TeV for (5 > 5.4. 

3.2. The SSC Spectrum 

We next discuss the SSC spectrum from the blob. 
We assume the accelerated electrons in the blob to 

have a power-law energy distribution with an index 
of p and the maximum Lorentz factor of 7m. Wc 
also assume that the magnetic field strength in the 
blob is B'. Since the synchrotron peak of the 
(p{ej) flux corresponds to the emission from electrons 
with the maximum Lorentz factor, we obtain 7m = 
2.93 X 106(VlO)~i/2(B'/o.lG)-i/2(gpj^/ioOkeV)i/2. Be- 
cause of 7Mebr ^ Smec"^, relativistic KN correc- 
tions to the Compton spectrum are very important. 
We estimate the lower limit to the observed gamma 
energy at which the KN effects should be consid- 
ered, ' Skn - ((5TOec2)Vebr = 13((5/10)2(ebr/2keV)-i 
GeV. In addition, the Gompton peak energy in the 
KN range is approximated by -B^^pk ~ S^mtt^cC^ = 
1.5((5/10)i/2(570.1G)-i/2(epi^/100keV)i/2TeV. Thus, 
the SSC photon spectrum oc E^^^~^^^^'^ for E^ < Ekn and 
cx E~P for -Bkn < E^f < -E^.pk- Comparing this spectrum 
with equation (1), wc find p = Fi ~ 1.6 — 1.7, which is 
not only consistent with the value of p shown by Kino et 
al. (2002) but is also in agreement with the index of a2 
obtained from the BeppoSAX observations by Plan et al. 
(1998). 

Figure 1 also presents the SSC spectra, whose high- 
energy segments are assumed to be the intrinsic spectra 
derived by DS. It can be seen from this figure that the 
flux of the externally scattered photons is higher than that 
of the SSC photons in the blob for an IGMF ^ 10"!^ G. 
Furthermore, the externally scattered photon flux is con- 
sistent with the existing observations by EGRET and large 
enough to be detected by GLAST. 

4. DISCUSSION AND CONCLUSIONS 

Using the intrinsic 7-ray spectrum of Mrk 501 during 
its 1997 flaring activity derived by DS as an example, we 
have predicted a new GeV emission component for this 
and similar HBLs, if their radiation reaches us through 
regions with a low magnetic field as expected in under- 
dense voids. This GeV emission is due to the IC scatter- 
ing of CMB photons by the electron/positron pairs pro- 
duced in interactions of high-energy photons with the cos- 
mic infrared-UV background photons. We summarize our 
findings: First, the GeV flux of the scattered CMB pho- 
tons is higher than that from the SSC process in the blob 
for IGMFs ^ IQ-^^G. For an IGMF ^ lO^^o Q, there is 
an observable spectral turnover, whose position is strongly 
dependent on the field strength. This spectral turnover is 
due to the fact that the lower energy emission would have 
a longer delay. Second, the typical duration of the GeV 
emission is always given by the variability timescale of the 
TeV gamma-ray flux for photon energes larger than E'tum, 

— 5/2 

and it is given by Ate oc E„. for photon energies smaller 
than i?turn- Third, the GeV emission predicted here is con- 
sistent with existing EGRET observations, and is strong 
enough to be detected by GLAST, for IGMFs ^ 10"!^ G 
and flare times ~ 1 day. Fourth, we have derived a generic 
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formula for the optical depth due to pair production in 
a relativistic blob, r^" oc 5^^"^"', where 5 is the Doppler 
factor and a is the index of the softer photon energy spec- 
trum. This optical depth is not only consistent with the 
previously used formula in the a = 1 case but can also be 
generalized to broader cases. 

We note that even if a source such as Mrk 501 is not 
in a void, primary TeV photons emitted from the source 
have such a long mean free path that electron/positron 
pairs may be produced in void regions. Thus, detections 
on the GeV emission would provide a sensitive probe for 
weak IGMFs, of consequence for cosmogonical and galactic 
dynamo theories. 

Our results are also relevant for the quiescent state of 
HBLs such as Mrk 501. Aharonian, Coppi, & Volk (1994) 
and Coppi & Aharonian (1997) discussed a similar model 
for very high energy emission (> 100 GeV) from blazars 
in the quiescent state. We here use the recent models for 
the intergalactic infrared- UV background radiation to dis- 
cuss lower energy emission from TcV blazars. According 
to Catanese et al. (1997), the Whipple telescope detected 
a flux of(l){E^ c:^ 1 TeV) = 8^3 x lO^^^ ergscm-^ 3-^ in the 
quiescent state of this HBL. The external pair-production 



optical depth r°^(£;^ = 1 TeV) ~ 0.5 and 0.6, and thus 
the flux of the externally scattered photons with an en- 
ergy of ~ 0.6 GeV is estimated as (j^iEj ~ 0.6 GeV) ~ 
(j){E^ ~ ITeV) X [exp(T:^^^) - 1] ~ 5 x lO'^^ gj-gscm-^ 
and 7 x 10^^^ ergs cm^^ s^^ for the "baseline" and "fast- 
evolution" models of DS, respectively. These values of 
the GeV emission flux are much larger than ~ 1.0 x 
10~^^ ergscm"^ s~^, the sensitivity of GLAST &i energy 
~ 0.6 GeV for steady sources in a one-year survey. There- 
fore, even in the quiescent state of Mrk 501, its externally 
scattered emission may also be detected by GLAST. It 
should be pointed out that such detections are possible 
for a magnetic deflection angle that is less than the blazar 
jet opening angle. This also requires that the field strength 
be below ^ 10^^^ G, similar to the flare case, for an open- 
ing angle ~ 0.1. Above this value of the magnetic field, 
the steady GeV emission may be suppressed. 
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Fig. 1. — High-energy 7-ray spectra of Mrk 501 for a flare with tvar = 0.5 days and rj = I. The spectral data are those measured by 
HEGRA {circles, Aharonian et al. 1999, 2001) and Whipple {squares, Krennrich et al. 1999), and their fitting is shown by line A. The 
EGRET upper limit at 100 MeV is from Catanese et al. (1997). Thick solid (B) ajid dashed (B') lines are synchrotron self-Compton intrinsic 
spectra, derived from the observed spectrum A after correction for absorption on the IR background "fast-evolution" and "baseline" models 
of DS, respectively. Thin lines (C and C') are the secondary photon spectra calculated in the text from the resulting pairs interacting with 
the CMB. The thin lines labelled with numbers — 5 correspond to intergalactic magnetic fields of zero, 10-20, 10-19, 10-18, 10-", and 
10-1® G, respectively. The thick dot-dash line represents the GLAST sensitivity computed for an exposure time of 0.5 days. 
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Fig. 2. — Same as Fig.l but for different flcire durations and "fast-evolution" IR background. The predicted signals are thin lines while 
three thick lines are the GLAST sensitivity, both for a given integration time. The solid, dot-dash and dashed lines correspond to tvar = 3 
days, 0.5 days and 1 hour for Big = 10~^* G. The thin dotted line is the signal for tvar = 3 days and Big = 10~^® G. 



